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Abstract. Wc study here the water-waves problem for uneven bottoms in a 
highly nonlinear regime where the small amplitude assumption of the Kortcweg- 
de Vries (KdV) equation is enforced. It is known, that for such regimes, a 
generalization of the KdV equation (somehow linked to the Camassa-Holm 
equation) can be derived and justified [7] when the bottom is flat. We gener- 
alize here this result with a new class of equations taking into account variable 
bottom topographies. Of course, the many variable depth KdV equations ex- 
isting in the literature are recovered as particular cases. Various regimes for 
the topography regimes are investigated and we prove consistency of these 
models, as well as a full justification for some of them. We also study the 
problem of wave breaking for our new variable depth and highly nonlinear 
generalizations of the KDV equations. 



1. Introduction 

1.1. General Setting. This paper deals with the water-waves problem for uneven 
bottoms, which consists in studying the motion of the free surface and the evolu- 
tion of the velocity field of a layer of fluid under the following assumptions: the 
fluid is ideal, incompressible, irrotationnal, and under the only influence of gravity. 
Earlier works have set a good theoretical background for this problem. Its well- 
posedness has been discussed among others by Nalimov [26], Yasihara [36], Craig 
[8], Wu [33], [34] and Lannes [22]. Nevertheless, the solutions of these equations 
are very difficult to describe, because of the complexity of these equations. At this 
point, a classical method is to choose an asymptotic regime, in which we look for 
approximate models and hence for approximate solutions. More recently Alvarez- 
Samanicgo and Lannes [2] rigorously justified the relevance of the main asymp- 
totical models used in coastal oceanography, including: shallow-water equations, 
Boussincsq systems, Kadomtscv-Petviashvili (KP) approximation, Green-Naghdi 
equations (GN), Serre approximation, full-dispersion model and deep-water equa- 
tions. Some of these models capture the existence of solitary water-waves and the 
associated phenomenon of soliton manifestation [17j . The most prominent example 
is the Korteweg-de Vries (KdV) equation [21] , the only member of the wider family 
of BBM-type equations |4] that is integrable and relevant for the phenomenon of 
soliton manifestation. The KDV approximation originally derived over flat bottoms 
has been rigorously justified in [5] [3D] H3 HE]- When the bottom is not flat, vari- 
ous generalizations of the KDV equations with non constant coefficients have been 
proposed EH] [32] E] E3 US EH ES El US E2] ■ One of the aims of this article is to 
justify the derivation of this Korteweg-de Vries equation with topography (called 
KDV-top). Another development of models for water-waves was initiated in order 

i 



2 



SAMER ISRAWI 



to gain insight into wave breaking, one of the most fundamental aspects of water- 
waves [12] . In 2008 Constantin and Lannes [7] rigorously justified the relevance of 
more nonlinear generalization of the KDV equations (linked to the Camassa-Holm 
equation [6] and the Degasperis-Proccsi equations [10] ) as models for the propoga- 
tion of shallow water-waves. They proved that these equations can be used to 
furnish approximations to the governing equations for water-waves, and in their 
investigation they put earlier (formal) asymptotic procedures due to Johnson [18] 
on a firm and mathematically rigorous basis. However, all these results hold for 
flat bottoms only. The main goal of this article is to investigate the same scaling 
as in [7] and to include topographical effects. To this end, we derive a new variable 
coefficients class of equations which takes into account these effects and generalizes 
the CH like equations of Constantin-Lannes [7]. The presence of the topography 
terms induce secular growth effects which do not always allow a full justification 
of the model. We however give some consistency results for all the models derived 
here, and then show that under some additional assumptions on the topography 
variations, the secular terms can be controlcd and a full justification given. 

1.2. Presentation of the results. Parameterizing the free surface by z = C,{t,x) 
(with x E R) and the bottom by z = —ho + b(x) (with ho > constant), one 
can use the incomprcssibility and irrotationality conditions to write the water- 
waves equations under Bernoulli's formulation, in terms of a velocity potential ip 
associated to the flow, and where ip(t, .) is defined on £l t = {(x, z), —ho + b(x) < 
z < £(t, x)} (i.e. the velocity field is given by v = V Xi2 </s) : 



where g is the gravitational acceleration, d n ip is the outward normal derivative at 
the boundary of the fluid domain. The qualitative study of the water-waves equa- 
tions is made easier by the introduction of dimcnsionlcss variables and unknowns. 
This requires the introduction of various orders of magnitude linked to the physical 
regime under consideration. More precisely, let us introduce the following quanti- 
ties: a is the order of amplitude of the waves; A is the wave-length of the waves; 
bo is the order of amplitude of the variations of the bottom topography; A/a is the 
wave-length of the bottom variations; ho is the reference depth. We also introduce 
the following dimcnsionlcss parameters: 



(1) 



' 3> + ^V = 0, in 
d n ip = 0, at 
' d t ( + d x (d x (p = d z ip, at 



z = C, 



z 



h + b, 



K dtv+\((d x p) 2 + (d z p) 2 ) + g( = at 




the parameter e is often called nonlinearity parameter; while /i is the shallowness 
parameter. We now perform the classical shallow water non-dimensionalization 
using the following relations: 
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so, the equations of motion ([I]) then become (after dropping the primes for the sake 
of clarity): 



(3) 



(id*ip + d*cp = 0, at -1 + /3b (a '> < z < e(, 

d z tp- n@ad x bWd x (p = at z = -l + a \ 

dtC ~ -(fied x (d x ip + d z (p) = 0, at z = e(, 

dtV+Ue(d x tp) 2 + -(d zV ) 2 ) + ( = Q at z = e£ 



where b^ a \x) = b(ax). Making assumptions on the respective size of e, (3, a, and 
\x one is led to derive (simpler) asymptotic models from ([3]). In the shallow- water 
scaling (fj. <C 1), one can derive (when no smallness assumption is made on s, (3 and 
a) the so-called Green-Naghdi equations (see [TBI [24] for a derivation and [2] for a 
rigorous justification). For one dimensional surfaces and over uneven bottoms these 
equations couple the free surface elevation ( to the vertically averaged horizontal 
component of the velocity, 
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(4) u[t,x) = — — — —j-r / d x <p(t,x,z)dz 

and can be written as: 

d t ( + d x {hu) = 0, 

J (l + ^T[h,0b^})dtU + d x ( + £ud x U 

+ iie{ - ^d x (h 3 (ud 2 x u) - {d x u) 2 ) + Z[h,l3b^}u} = 
where h = l + e(- a) and 

T[h,f3b^]W = -^d x (h 3 d x W) + ^d x (h 2 d x b^)W + f3 2 h(d x b^) 2 W, 
while the purely topographical term $s[h, f3b^]u is defined as: 

= ^[d x {h 2 {ud x ) 2 b^)-h 2 {{ud 2 x u)-{d x u) 2 )d x b^)\ 

+(3 2 {{ud x ) 2 b^)d x b^ a \ 

If we make the additional assumption that e <C 1, (3 <C 1 then the above system 
reduces at first order to a wave equation of speed ±1 and any perturbation of the 
surface splits up into two components moving in opposite directions. A natural issue 
is therefore to describe more accurately the motion of these two "unidirectional" 
waves. In the so called long- wave regime 

(6) n « 1, e = 0(/i), 

and for flat bottoms, Korteweg and de Vries [3JJ found that say, the right-going 
wave should satisfy the KDV equation: 

(7) u t + u x + \euu x + ^u xxx = 0, 

2 o 

and (C = u + 0(e, n)). 

At leading order, this equation reduces to the expected transport equation at speed 
1. It has been noticed by Benjamin, Bona, Mahoney [I] that the KDV equation 
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belongs to a wider class of equations. For instance, the BBM equation first used 
by Peregrine [27j . and sometimes also called the regularized long- wave equation, 
provides an approximation of the exact water-waves equations of the same accuracy 
as the KDV equation and can be written under the form: 

3 1 

(8) u t + u x + -euu x + /j,(Au xxx + Bu xxt ) = with A — B = -. 

I o 

For higher values of e, the nonlinear effects are stronger; in the regime 

(9) e = 0(v^)> 

the BBM equations ([8]) should be replaced by the following family (see [71 ITS"]): 
3 

(10) u t + u x + -euu x + fi(Au xxx + Bu xxt ) = E[i(Euu xxx + Fu x U xx ^) 

(with some conditions on A, B, E, and F) in order to keep the same 0{y?) accuracy 
of the approximation. However, all these results only hold for flat bottoms; for the 
situation of an uneven bottom, various generalizations of the KDV equations with 
non constant coefficients have been proposed [20l l32l [TT | [25] . [TB I [35] . [25 1 fl4 l H^ l [29] ■ 
We justify in this paper the derivation of the generalized KDV equation and also 
we show that the correct generalization of the equation (|10jl under the scaling ((9]) 
and with the following conditions on the topographical variations: 

(11) /3a = 0(jx), pa 2 =0{fi 2 ) /3as = 0{fj?), 
is given by: 

u t + cu x + \c x u + \euu x + fi(Au xxx + Bu xxt ) 

(12) = euEuu xxx + en(d x {^u)u xx + u x d£(^uf) 

where c = y/l — [3b^ and A, E, F differ from the coefficients A, E, F in (fTU)) 
because of topographic effects: 

A = Ac 5 - Be 5 + Be 

E = Ec 4 --Bc 4 + -B 

2 2 

q Q 
F = Fc 4 --Bc 4 + -B. 

2 2 

Notice that for an equation of the family (fl2~|) to be linearly well-posed it is necessary 
that B < 0. In Sec. we derive asymptotical approximations of the Green-Naghdi 
equations over non flat bottoms: equations on the velocity are given in Sect. 12.11 
and equations on the surface elevation are obtained in Sect. 12. 2\ for these equations, 
inconsistency results are given (see Definition [1]) . In Sect. 12.31 the same kind of 
result is given in the (more restrictive) KdV scaling in order to recover the many 
variable depth KdV equations formally derived by oceanographers. Section [3] is 
devoted to the study of the well posedness of the equations derived in Section [5] 
Two different approaches are used, depending on the coefficient B in Q12p : £13.11 
deals with the case B < (in that case, further investigation on the breaking of 
waves can be performed, see fc|3.2p and ij3.3l treats the case B = 0. While secular 
growth effects prevent is from proving i? s -consistency (see Definition [2]) for the 
models derived in Section [21 we show in Section [4] that such results hold if one 
makes stronger assumptions on the parameters. A full justification of the models 
can then be given (see Th. [4|. 
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2. Unidirectional limit of the Green-Naghdi equations over uneven 
bottom in the ch and kdv scalings 

We derive here asymptotical approximations of the Green-Naghdi equations over 
non flat bottoms in the scalings (fTTjl and ([9]). We remark that the Green-Naghdi 
equations can then be simplified into 
(denoting h = 1 + e( - /3b a ): 

( Q + [hu} x = 

(13) < . ll o . 2 \i 

[ ut + Cx + £uu x = — [h [u xt + euu xx - eu x )\ x , 

where 0(li 2 ) terms have been discarded. 

We consider here parameters e, j3, a and fx linked by the relations 

(14) e = O(Vm), Pet = 0(e), f3a = 0{n), (3a 2 = 0{ll 2 ), (3ae = 0(li 2 ) 

(note that in the case of flat bottoms, one can take (3 = 0, so that this set of 
relations reduce to e = O(^fji)). 

Equations for the velocity u are first derived in §2.1 and equations for the surface 
elevation £ are obtained in §2.2. The considerations we make on the derivation of 
these equations are related to the approach initiated by Constantin and Lannes [7] . 
In addition, in §2.3 we recover and justify the KDV equation over a slowly varying 
depth (formally derived in [501 1521 ITT] ) . 

2.1. Equations on the velocity. If we want to find an approximation at order 
0(fj, 2 ) of the GN equations under the scalings (fT4")) , it is natural to look for u as a 
solution of (fT2")) with variable coefficients A, B, E, F to be determined. We prove 
in this section that one can associate to the solution of (|12[) a family of approximate 
solutions consistent with the Green-Naghdi equations (fT"3|) in the following sense: 

Definition 1. Let p be a family of parameters 9 = (e,/3,a, ll) satisfying |1^[ ). A 
family (Q® ,u ( - e )e 6p is L°° -consistent on [0, with the GN equations U3\) . if for all 
9 e p (and denoting h 6 = 1 + eC, e - a) ), 

£ (4) 2 )]x + A 2 e 

with {r{,rl)e ep bounded in L°°([0, ~] x K). 

Remark 1. The notion of L°° -consistency is weaker then the notion of H s -consistency 
given in §4 (Definition^) and does not allow a full justification of the asymptotic 
models. Since secular growth effects do not allow in general an H s -consistency, we 
state here an L°° -consistency result under very general assumptions on the topog- 
raphy parameters a and (3. H s -consistency and full justification of the models will 
then be achieved under additional assumptions in §^. 

The following proposition shows that there is a one parameter family of equations 
(fT2"]l L°°-consistent with the GN equations (fT5]) . (For the sake of simplicity, here and 
throughout the rest of this paper, we take an infinity smooth bottom parameterized 
by the function b). 



Cf + [h e u% = 
u\ + C + eu e u 



„ 2 r e 
P '1 



3h e 



{(hy(u xt + eu» 
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Proposition 1. Let b 6 H°°(M.) and p £ R. Assume that 

1 3 1 e. 9 23 
A = p, B = p h = — p , h = — p . 

y ' F 6 2 F 6' 2 F 24 

TTien: 

• For all family p of parameters satisfying Ji^[ ). 

• For all s > /arge enough and T > 0, 

• For a// bounded family {u e )ee P € C([0, f ], H" S (R)) solving IW) . 
the familly (<^ e ,u e )g ep with (omitting the index 9) 



/-ir-\ t" ^ f ^2M4 4r^ ^ 2 1 

15 C := cu H — / c T u H — it H — c u T + — sue \—uu TX -\ uJ, 

zs L°° -consistent on [0, -j-] with the GN equations U3\). 

Remark 2. If we take b — -i.e i/ioe consider a fiat bottom-, then one can recover 
the equation (7) of [7] and the equations (26a) and (26b) of [18j with p = — 
p = g respectively. 
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Proof. For the sake of simplicity, we denote by 0(/x) any family of functions (/ )eea 

such that —f 6 remains bounded in L°°([0, ~], H r (R)) for all 6 G p, (and for possibly 
t 1 

different values of r). The same notation is also used for real numbers, e.g e = O(p), 
but this should not yield any confusion. We use the notation Oi=« (p) if — f e remains 

bounded in L°°([0, ^] x R). Of course, similar notations are used for 0(fi 2 ) etc. To 
alleviate the text, we also omit the index 9 and write u instead of u . 
Step 1. We begin the proof by the following Lemma where a new class of equation 
is deduced from (fT2|) . The coefficients A, B, E, F in this new class of equations 
are constants (as opposed to A, E and F in ([12"]) that are functions of x). 

Lemma 1. Under the assumptions of Proposition [7J there is a family (R e )e^ p 
bounded in L°°([0, -^],iif r (R)) (for some r < s) such that (omitting the index 9) 

(16) u t + cu x + \c x u + \euu x + fic 5 Au xxx + iiBd x (c A u xt ) 

= e/j,c i Euu xxx + eh\F{c a u) x u xx + sfj,^Fu x (c i u) xx + p?R. 

Proof. Remark that the relation a(3 = 0(fj,) and the definitions of A, B, and E in 
terms of A, B, E and F imply that 

3 3 

H(-Bc 5 + Bc)u xxx + efi(--Bc 4 + -B)uu xxx 

+ei*dx { {— 2 — — J uj u xx + epiu x d 2 x ^-^ — a_j u j 

3 

= -fJiBd x (c(c 4 - l)d x u x ) - -^eBd x ((c A - l)d x (uu x )) 



3 1 

(c 4 - l)d x (cu x + -euu x ) + Oip 2 ) 



= -fiBd x 
= fiBd x {{c 4 - l)u xt ) + Q{p 2 ) 
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the last line being a consequence of the identity ut = —{cu x + ^euu x ) + 0(/i) 
provided by (|12p since we have 

\c x u\ H r = | - —Pa$?>v\ HT < Cst aj3\d x b^\ wlr]+1:00 \u\ H r = 0(/3a) = O(fi) 

where [r] is the largest integer smaller or equal to r. The equation (|12p can thus 
be written under the form: 

u t + cu x + \c x u + \euu x + /_ic 5 Au xxx + ^Bd x (c A u xt ) 
= £\ic A Euu xxx + ei±\F(c 4 u) x u xx + £^Fu x (c i u) xx + 0(^ 2 ), 

which is exactly the result stated in the Lemma. □ 
If u solves (fl2ll one also has 



3 3 

(17) ut + cu x + -euu x = --c x u + 0(fi) 

= O00, 

Differentiating (|17p twice with respect to x, and using again the fact that c x f = 
0((3a) = O(n) for all / smooth enough, one gets 

It is then easy to deduce that 

3 

c 5 u xxx = ~c 4 u xxt - -c A ed 2 x (uu x ) + 0(fi) 

3 

= -d x (c 4 u xt ) - -£c A {uu xxx + 3u x ii xx ) + O(p) 

so that we can replace the c 5 u xxx term of (|16p by this expression. By using Lemma 
[TJ one gets therefore the following equation where the linear term in u xxx has been 
removed: 

3 3 

(18) u t + cu x + -c x u + -euu x + [ic A au xxt = £fxc [euu xx + du 2 x ] x + 0{^ 2 ) 
with a = B-A,e = E+\A,d = \{F + 3A- E). 

Step 2. We seek v such that if C := cu + £v and u solves (TT2]) then the second 
equation of (fT3"|) is satisfied up to a 0(/i 2 ) term. This is equivalent to checking that 

Ut + [cu + £v] x + £uu x = fi(l - (3b < - a) )£d x (cu)u xt + -((1 - [3b {a) ) + £cufu xxt 

+ £ JL(l _ p b ^) 2 (uu xx - u 2 x ) x + 0(fi 2 ) 



LL ( 2 C 4 \ 

— c u xxt + £n(c 3 u x u xt + -c 3 uu xxt + —{uu xx - u x ) x J 
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where we used the relations 0(e 2 ) = O(n), 0((3a) = 0(/i) and the fact that 
c 2 = 1 — . This condition can be recast under the form 

ev x + [u t + cu x + \c x u + \euu x + ^au xxt - e^[euu xx + du%] 

IS 4 .1. 

= -^c x u + -uu x + fic (a + -)u xxt 

+e^c 3 (u x u xt + ^uu xxt + c[(- - e)uu xx - (- + d)u 2 x ] x ^ + 0{^ 2 ). 

Since morever one gets from (fT7j) that u x t = —cu xx + 0(/j,,e) and u xx t = —cu xxx + 
0(iJ,,e), one gets readily 

ev x + [u t + cu x + \c x u + \euu x + ^au xxt - e»[euu xx + du x ] x ] 

1 £ 1 1 1 

= 2 C z u + 2 UUx + /^ c4 ( a + g) - "^* - e/ic 4 [(e + + (d + -)u x } x + 0{ji 2 ). 

From Step 1 , we know that the term between brackets in the lhs of this equation is 
of order 0(/i 2 ) so that the second equation of (flB"]) is satisfied up to 0(/j, 2 ) terms if 

(19) ev x = ic^u+|wu x + ^c 4 (a+i)M XKt -e^c 4 [(e + i)uu^ + (d+i)u 2 .] x + 0(^ 2 ). 
At this point we need also the following lemma 

Lemma 2. With u and b as in the statement of Proposition [IJ the mapping 
(t, x) — ► /_ c x u dx is well defined on [0, x R. 
Morever one has that 



c x u dx 



L~([0,Z]xR) 



Proof. We used here the Cauchy-Schwarz inequality and the definition c 2 = 1— 
to get 

fX 

c x u dx < Cst af3\{b x Y a) \ 2 \u\2 < Cst Vot0\b x \ 2 \u\ 2 < oo. 

3 

It is then easy to conclude the proof of the lemma. □ 

Thanks to this lemma there is a solution v G C([0, — ] x K) to ([IT))) , namely 

1 e 1 1 1 

(20) ev = - J c x u + -u 2 + fic 4 (a+ -)u xt - epic 4 [(e + -)uu xx + (d+ -)u 2 x ]. 

Step 3. We show here that it is possible to choose the coefficients A, B, E, F 
such that the first equation of (fT^|) is also satisfied up to Ol^{^l 2 ) terms. This is 
equivalent to checking that 

(21) [cu + ev] t + [(1 + s{cu + ev) - (3b {a) )u] x = 

Remarking that the relations 0((3a) = 0(/i), 0{(3a 2 ) = 0(fi 2 ), and 0((3ae) = 
0((j?) imply that 

I r x i 

■Z I C x U t = -~CC X U + L <*=((J, ), 
^ J-oo ^ 
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one infers from (l20l) that 



ev t = c x u t + ^uu t + fic 4 {a + ^)u xtt - efic 4 [(e + ^)uu xx + (d+^)ul} t 

= -^cc x u - |u(cu a + y UUx + M^ouxast) - ,uc 4 (a + ^)<9 2 t (cu x + 
+£/^c 5 [(e + -)««„. + (d + + 0(^ 2 ) + L =c (/?) 

1 1 2 3 2 , Is 

= --cc x u - e-cuu x - e -u u x - n{a + - )c u xxt 
5 5 

+e/ic 5 [(2a + e+ -)uu xx + (-a + d + l)u 2 x ] x + 0(/i 2 ) + Ol°°(// 2 ). 
4 

Similarly, one gets 

3 1 

e 2 [vu] x = e 2 -u 2 u x - e^(a + -)[uu xx ] x + 0{ii 2 ) + L ^(fj, 2 ), 

so (|21jl is equivalent to 

cu t + evt + c 2 u x + 2cc x u + 2ecuu x + e 2 [vu] x = Ol°° (a* 2 )- 

Multiplying by -, we get 

c 

u t + cu x + ^c x u + e^uu x - ^c 4 (a + -)u xxt 
1 5 

= e^c 4 [-(e + a + -)uu xx - (-a + d + l)w 2 ]:r + O l ^(^ 2 )- 

Equating the coefficients of this equation with those of (fT5|) shows that the first 
equation of (fl~3"|) is also satisfied at order Ol°°{^ 2 ) if the following relations hold: 

1 1 , 19 

a = , e — , a = , 

6' 6' 48' 

and the conditions given in the statement of the proposition on A, B, E, and F 
follows from the expressions of a, e and d given after (fl~8"|) . □ 

2.2. Equations on the surface elevation. Proceeding exactly as in the proof of 
Proposition [TJ one can prove that the family of equations on the surface elevation 



(22) +M(iC^ + BC rat ) = epECCxxx +ev(d x {-0( xx + C^ 2 (f C)), 

where 



I = ,4c 5 - Sc 5 + Be 

E = Ec 3 - -Be 3 + —B 
2 2c 

F = Fc 3 - -Be 3 + —B, 
2 2c ' 

can be used to construct an approximate solution consistent with the Green-Naghdi 
equations: 
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Proposition 2. Let b 6 H°°(M.) and ggi. Assume that 



A = q, B = q-- E = --q--, F = --q-—. 
*' 1 6 2 y 6 2 y 24 



Then: 



• For all family p of parameters satisfying {1$ , 

• For all s > large enough and T > 0, 



i<br aZZ bounded family (( e )g ep £ C([0, — ],ff s (M)) solving 



the familly (<^ e ,u e )g ep with (omitting the index 0) 



- M Ic 3 C, t + £MC 2 [icCxx + iCx])) 



6 ' ' (i 48 

is L°° -consistent on [0, iot£/i £/ie GiV equations Ul 

Remark 3. If we take 6 = -i.e if we consider a fiat bottom-, then one can recover 
the equation (18) o/[7]. 

Remark 4. Choosing q = A, a = e and (3 = /i 3 / 2 the equation k22\) reads after 
neglecting the 0(/i 2 ) terms: 

Ct + cCx + \c x C + |eCC« - ^ 2 c 2 a + |^ 3 C 3 C* 

(24) "I" (Casasas Cxat) = TJ^^A'CCCxaja; 2^^^), 



it is more advantageous to use this equation \2J$ to study the pattern of wave- 
breaking for the variable bottom CH equation (see $3.2 below). 

Proof. As in Proposition [TJ we make the proof in 3 steps (we just sketch the proof 
here since it is similar to the proof of Proposition [TJ). 
Step 1. We prove that if C, solves (f22j) one gets 



(25) Ct + <x + \c x u + ±eC( x - + T^e 3 C 3 a + tic 4 a( xxt 

= e^KCra + dQ] x + 0(/i 2 ) 
with a = B - A, e = E + ^A, d= \{F + 3A - E). 

Step 2. We seek v such that if u := — (C + ev) an d C solves (f2"2")l then the first 

c 

equation of (Tl3|) is satisfied up to a Ol^(^ 2 ) term. Proceeding as in the proof of 
Proposition [TJ one can check that a good choice for v is 

+fj,c 3 a( xt - £^c 2 [eQC, xx + dQl\. 

Step 3. We show here that it is possible to choose the coefficients A, B, E, F such 
that the second equation of (fl"3"|) is also satisfied up to Ol°°{h 2 ) terms. Replacing 

u by — (C + £ y ) with v given by (|26[) . one can check that that this condition is 
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equivalent to 

1 3 3 3 1 

Ct + < x + t;C x ( + —eCCx - ^ C Cx + —^e 3 C 3 Cx - ^{a + -)( xxt 

2 2c 8c J 16c 3 

= £/ic 3 [-(e + a- ijCCx* -{\a + d+ ^)(% + L ~ ( M 2 ). 

Equating the coefficients of this equation with those of (|25p shows that the second 
equation of ((T3|) is also satisfied at order Ol°°(a* 2 ) if the following relations hold: 

1 1 i 1 

a = — , e = — , a — , 

6' 6' 48' 

and the conditions given in the statement of the proposition on A, B, E, and F 
follows from the expressions of a, e and d given after (|2"5|) . □ 

2.3. Derivation of the KdV equation in the long-wave scaling. In this sub- 
section, attention is given to the regime of slow variations of the bottom topography 
under the long- wave scaling e = 0(/i). We give here a rigorous justification in the 
meaning of consistency of the variable-depth extensions of the KdV equation (called 
KdV-top) originally derived in [301 [331 [TT] . We consider the values of e, /?, a and 
fj. satisfying: 

(27) e = C>(/i), a/3 = 0(e), a 2 /3 = 0(e 2 ). 

Remark 5. Any family of parameters 9 — (e,/3, a, /i) satisfying also satisfies 

m ■ 

Neglecting the 0(/z 2 ) terms, one obtains from (|13[) the following Boussinesq 
system: 

f Ct + [hu] x = 

(28) \ /i 4 

[ u t + Q x + euu x = —c u xxt , 

where we recall that h = 1 + e( — fib^ and c 2 = 1 — [3ty- a >. The next proposition 
proves that the KdV-top equation 

3 1 1 

(29) Ct + < x + — eCC* + ^c b C, xxx + -c x ( = 0, 

2c D 2 



is F°°-consistent with the equations 

Proposition 3. Let b G H°°(R). Then: 

• For all family p' of parameters satisfying ^P, 

• For all s > Zarge enough and T > 0, 

• For a/Z bounded family (( e )g ep > G C([0, j], if s (R)) solving 

the familly {(f ,u e )e£ P > with (omitting the index 6) 



is L°° -consistent on [0, u;ii/j i/*e equations A28\) . 

Remark 6. Similarly, one can prove that a family (C e , u e )g^ p > with u e solution of 
the KDV -top equation 

3 1 3 

(30) u t + cu x + -euu x + -nc 5 u xxx + -c x u = 0, 

zo 2 
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and (f given by 
(31) 



C:=cu+- 



£ 2 M <i 

c T u H — u c Ur, 

4 6 



is L°° -consistent with the equations i28\ ). 

Proof. We saw in the previous subsection that if ((^Wp is a family of solutions of 
(|22|) . then the family (£ e , u e )e ep with is given by (|23|) is L°°-consistent with the 
equations f) 1 3|) . Since 6 = (e, /3, a, fi) G p' C p then by taking q = i, we remark that 
the equations (|^|) and (f!2l?)) are equivalent in the meaning of L°°-consistency and 
the systems (flU)) and (f^5|) arc also, so it is clear that (C , u e )g^ P ' is L°°-consistent 
with the equations (1281) . □ 



3. Mathematical analysis of the variable bottom models 

3.1. Well-posedness for the variable bottom CH equation. We prove here 
the well poscdness of the general class of equations 

(1 - /.imdl)u t + cu x + kc x u + ^ e j f 3 u 3 u x 



(32) 



6fi 



2 



+ d x (h 2 u)u xx + u x d x {h 2 u) 



where m>0,/ceIR,Jisa finite subset of N* and fj = fj(c), g = g(c), hi = hi(c) 
and h 2 = h 2 (c) are smooth functions of c. We also recall that c = y/l — /3b( a \ 
Example 1. Taking 

3 

~ 2' 

3 



m = —B, 



2' 



g(c) = Ac 5 - Be 5 + Be, hx{c) = c 4 E - -Be 4 



3 -B, 
2 



/12(c) 



Fc 4 - \Bc 4 + \B 



the equation \32]) coincides with \12] 
Example 2. Taking 



m = —B, 



/i(c) = 



2c' 



2' 

/ 2 (c) = - 



J = {1,2, 3}, 
3 

8c3' 



/3(c) 



16c 5 ' 



g(c) = Ac 5 - Be 5 + Be, 

Fc 3 - §5c 3 + l c B 
/12(c) = 2 — —, 



h^c) = c 3 E--Bc 3 



2c 



-B, 



the equation h32]) coincides with 

More precisely, Theorem Q] below shows that one can solve the initial value 
problem 

(1 - fimd x )u t + cu x + kc x u + ^ e 3 fjU j u x + [igu xxx 
= S(i hiuu xxx + d x (h 2 u)u xx + u x dl(h 2 u) , 



(33) 
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on a time scale 0(1/ e), and under the condition m > 0. In order to state the result, 
we need to define the energy space X s (s <G R) as 

X S+1 (R) = H S+1 (R) endowed with the norm |/|| s+ i = \f\ 2 H . + fJ,m\d x f\ 2 Hs . 

Theorem 1. Let m > 0, s > | and b e H°°(M). Let also p be a family of parame- 
ters 9 = (s,f3,a,fj,) satisfying Then for all u° G H S+1 (M), there exists T > 
and a unique family of solutions (u e )g ep to Ii33\) bounded in C([0, — ];X S+1 (R)) n 

Proof. In this proof, we use the generic notation 

C = C(e,fi,a,f3,s, \b\ H °) 

for some a > s + 1/2 large enough. When the constant also depends on |u|x s +i, 
we write Cd^lx^+O- Note that the dependence on the parameters is assumed to 
be nondecreasing. 

For all v smooth enough, let us define the "linearized" operator C(v,d) as 

C(v,d) = (1 — fj,md 2 )d t + cd x + kc x + e^fjv^d x + [igd\ 
-e^hxvdl + (h 2 v) x dl + (h 2 v) xx d x ]; 

for the sake of simplicity we use the convention of summation over repated indexes, 
Sjej^/j^ = ^ 'fjW '■ To build a solution of (j3"2")) using an iterative scheme, we 
have to study the problem initial value 



(34) 



C(v,d)u = ef, 
u U=o = u °- 



If v is smooth enough, it is completely standard to check that for all s > 0, / G 
Lj oc (R+;H s (R x )) andu G H S (R), there exists a unique solution u G C"(M+; H s+1 (R)) 
to (f34|) (recall that m > 0). We thus take for granted the existence of a solution 
to ([34]) and establish some precise energy estimates on the solution. These energy 
estimates are given in terms of the | ■ \x s + 1 norm introduced above: 



\ 2 H s + fxm\d x u\ 2 H s 



Differentiating eAt |it| Xs+ i with respect to time, one gets using the equation 
and integrating by parts, 

^ xt d t (e-" xt \u\ 2 xs+1 ) = -yM| s+1 - (A s (cd x u), A s u) - k(A s (ud x c), A s u) 

+e(A s f, A s u) - e j (A s (f jv ] d x u) , A s u) - (,(A s (gd x u), A s u) 
-eft{A'(hivd£u),A'u) - e^(A s ((h 2 v) x d x u), A s d x u) 7 

where A = (1 — d 2 ) 1 / 2 . Since for all constant coefficient skewsymmctric differential 
polynomial P (that is, P* = — P), and all h smooth enough, one has 

(A s {hPu),A s u) = ([A s ,h]Pu,A s u)- -([P,h]A s u,A s u), 
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we deduce (applying this identity with P = d x and P = d x ), 
\e^d t {e-^\u\ 2 x^) = + £ ( A V, A s u) 

-([A s , c]d x u, A s u) + -((d x c)A s u, A s u) - k{A s (ud x c), A s u) 

'( A". J ,r r <) r u. \ ~ ti) + j(d x (fjV j )A s u,A s u) 
3 

-fj,([A s ,g]d 2 u- -g x A s d x u-g xx A s u,A s d x u) - n([A s ,g x ]dlu 7 A s u) 
3 

-en([A s ,h lV }dlu~- -{h lV ) x A s d x u- [hiv) xx A s u, A s d x u) 
(35) — e^([A s , (hiv) x ]d x u, A s u) — en(A s ((h 2 v) x d x u),A s d x u). 
Note that we also used the identities 

[A s , h]d 3 x u = d x ([A s , h]d 2 x u) - [A s , h x ]d 2 x u 

and 

-(h xxx A s u,A s u) = -(h xx A s u,A s u x ). 

The terms involving the velocity c (second line in the r.h.s of |35|) ) are controled 
using the following lemma: 



Lemma 3. Let s > 3/2 and (e,f3,a,fi) satisfy \14\) - Then there exists C > such 
that 

(36) \([A s ,c]d x u,A s u)\ < eC\u\ 2 Xs+1 

(37) \((d x c)A s u,A s u)\ < eC\u\\ s+1 

(38) \{A s {ud x c),A s u)\ < eC\u\ 2 Xs+1 . 

Proof. • Estimate of | ([A s , c]d x u, A s u) |. One could control this term by a standard 
commutator estimates in terms of |c K |^ s -i; however, one has Ic^l^s-i = 0(yfa(3) 
and not 0(a(3) as needed. We thus write 

([A s ,c\d x u, A s u) = (([A s ,c] - {A s ,c})d x u 7 A s u) + ({A s ,c}d x u,A s u), 

where for all function F, {A s , F} stands for the Poisson bracket, 

{A S ,F} = -sd x FA s ~ 2 d x . 

We can then use the following commutator estimate ([53], Theorem 5): for all F 
and U smooth enough, one has 

Vs > 3/2, |([A S ,F] - {A S ,F})U\ 2 < Cst |^F| H .|^._ 2 . 

Since \d 2 c\n= = 0(a(3) = 0(e), we deduce 

|(([A s ,c] - {A s !C })^,A s n)| < C/3a\u\ 2 X3+1 < Ce\u\% a+1 . 

Morever 

\({A s ,c}d x u,A s u)\ < \-sd(c 4 )A s - 2 d 2 M 2 \u\ Xs+1 < C\d(c 4 )Uu\ 2 XB+1 < eC\u\ 2 Xs+1 , 
and ((36|) thus follows easily. 

• Estimate of | ((d x c)A s u, A s uj |. Since \u\h s < \u\x*+i, one has by Cauchy-Schwarz 
inequality 

((d x c)A s u,A s u) < \d x c\ 
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therefore, by using the fact that |9a;c|oo < Cst af3\d x b\oc and af3 = 0(e), one gets 
easily 

((d x c)A s u,A s u) <eC\u\ 2 Xs+1 . 

• Estimate of | (A s (ud x c), A s u) |. By Cauchy-Schwarz we get 

\(A s (ud x c),A s u)\ < \ud x c\ H s\u\ H >. 

We have also that 

\ud x c\ H s < |9 x c| V i / [ s j + i,oo|u|x=+i 

where [s] is the largest integer smaller or equal to s (this estimate is obvious for 
s integer and is obtained by interpolation for non integer values of s.) Using this 
estimate, and the fact that (3a = 0(e), it is easy to deduce ([38| . □ 

For the terms involving the fj, (third line in the r.h.s of (|35|) ) we use the controls 
given by: 

Lemma 4. Under the assumptions of Theorem^ one has that 

(39) \e^([A s jy}d x u,A s u)\ < eC(\v\ x .+i)\u\ i x . +1 

(40) \e^(d x (f^)A s u,A s u)\ < eC(\v\ Xs+l )\u\ 2 Xs+1 . 

Proof. • Estimate of e 3 ([A s , fjV 3 ]d x u, A s u). By Cauchy-Schwarz we get 

\e 3 ([A s ,f jV 3 ]d x u,A s u)\ < e 3 \[A s ,f 3 v 3 ]d x u\ 2 \u\x^- 

Let us recall here the well-known Calderon-Coifman-Mcycr commutator estimate: 
for all F and U smooth enough, one has 

Vs > 3/2, \[A S ,F]U\ 2 < Cst\F\ H s\U\ H s-i; 

using this estimate, it is easy to check that one gets (|39l) . 

• Estimate of e 3 (d x (fjV 3 )A s u, A s u) . It is clear that 

\(d x (f 3 v 3 )A s u,A s u)\ < Id^f^looHx.+i- 
Therefore (|40|) follows from the continuous embedding H s C W l o ° (s > 3/2). □ 

Similarly, the terms involving g (fourth line in the r.h.s of (|35| ) are controled 
using the following lemma: 

Lemma 5. Under the assumptions of Theorem^ one has that 

(41) | - f i([A s , g]d 2 x u - ^g x A s d x u - g xx A s u, A s d x u) 



^([A s ,g x ]dlu,A s u) <eC\u\ 2 x 



2 

+ 1 ■ 



Proof. Since \u\hs < \u\ Xs +i and ^/Jl\d x u\H' < -^\u\ Xs +i , by using the Cauchy- 
Schwarz inequality and proceeding as for the proofs of Lemmas[3l|4]one gets directly 

(SB. □ 



Finally to control the terms involving hi, (fifth and sixth lines in the r.h.s of 
([35l) ) let us state the following lemma: 
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Lemma 6. Under the assumptions of Theorem^ one has that 

(42) | - ^([A s , h x v\dlu - ^(h lV ) x A s d x u - {h lV ) xx K s u, A s d x u) 

- ¥ ([A 8 ,(M I ]^,A 8 «)| <eC(\v\ x , + i)\u\ x , +1 , 

(43) \-efi(A s ((h 2 v) x d x u),A s d x u)\ < eC(\v\ x * + i)\u\ 2 X3+1 . 
Proof. We remark first that 

\n({hiv) xx A s u, A s d x u) \ < ^(VmSxCM))! 
since s — 1 > i, so by using the imbedding H S ^ 1 (M.) C L°°(K) we get 

|9x(v7*^(/iiu))|oo < C(|?;| X <,+i). 
Proceeding now as for the proof of Lemma [5] one gets directly (|42[) and (|43]) . □ 

Gathering the informations provided by the above lemmas, we get 

e £At a t (e- £At | U |x s+1 )<e(C(|«|x^0-A)| W |^ +1 +2 £ |/U s+1 | u | xs+1 . 

Taking A = At large enough (how large depending on sup tg j tj C(\v(t)\ x „+i) to 
have the first term of the right hand side negative for all t e [0, j], one deduces 

Vt e [0,|], 9 4 (e- £A -*| U || s+1 ) < 2Ee-^t\ f \ xs+l \ u \ xs+l . 
Integrating this differential inequality yields therefore 

V* G [0, -], \u(t)\ X s +1 < e eXTt \u°\ X s+i + 2e [ e eXT < t -''>KO|jf+*#'- 

Thanks to this energy estimate, one can conclude classically (see e.g. pQ) to the 
existence of 

T = T(\u°\ xs+l ) >0, 

and of a unique solution u € C([0, — ]; X s+1 (R d )) to (|33[) as a limit of the iterative 
scheme 



w = u , and \/n e N, 



£{u n ,d)u n+1 = 0, 



u] l+1 = u°. 

|t=o 



Since u solves (f32j) . we have £(u,d)u = and therefore 

(A s_1 (l - iimdl)d t u, A s ~ 1 d t u) = -(A s ~ 1 M(u, d)u, A s_1 <9 t u), 
with 9) = £(u, 9) — (1 — fimd x )dt- Proceeding as above, one gets 

la^l*. <c(\u°\ X s + i,\u\ X s+i), 

and it follows that the family of solution is also bounded in C 1 ([0, — ]; X s ). □ 
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3.2. Explosion condition for the variable bottom CH equation. As in the 

case of flat bottoms, it is possible to give some information on the blow-up pattern 
for the equation (|2"4")) for the free surface. 

Proposition 4. Let b 6 H°°(M.), Co <= H 3 (R). If the maximal existence time 
T m > of the solution of with initial profile C(0, •) = Co is finite, the solution 
C € C([0,T m );H 3 (R)) n C^j0,T m );H 2 (R) is such that 

(44) sup {|CM|} <oo 

te[o,T m ),x&s. 

and 

(45) sup{Cc(i, x)} t oo as 1 1 T m . 

Remark 7. It is worth remarking that even though topography effects are introduced 
in our equation \21$ , wave breaking remains of 'surging' type (i.e the slope grows 
to +oo ) as for flat bottoms. This shows that plunging breakers (i.e the slope decays 
to -oo ) occur for stronger topography variations then those considerd in this paper. 

Proof. By using the Theorem [T] given Co G H 3 (R), the maximal existence time of 
the solution C(i) to ([24]) with initial data £(0) = Co is finite if and only if |C(i)|ff 3 (K) 
blows-up in finite time. To complete the proof it is enough to show that 

(i) the solution C(t) given by Theorem [T] remains uniformly bounded as long as 
it is defined; 

and 

(ii) if we can find some M = M(Co) > with 

(46) £ x (t,x)<M, xeR, 

as long as the solution is defined, then |C(^)li/ 3 (R) remains bounded on [0, T m ) . 
Remarking that 

/ [«* + \c x QQ dx = - I [(cC)xC " \c x ( 2 } dx = - f [(cC + \c x QQ dx, 
we can deduce 

(47) f [(cCx + JcxOC] dx = 0. 

JR * 

We have also that 

(48) [[(CC*)Qdx = 0, VzeN*. 

JR 

Item (i) follows at once from gT]), (J48J) and the imbedding F 1 (R) C L°°(M) since 
multiplying (|2~4")) by C and integrating on K yields 

(49) d t ([[e + ^nUl\dx)=Q. 

K JR L * JR ' 

To prove item (ii), as in [7] we multiply the equation (|24[) by Cxxxx and after 
performing several integrations by parts we obtains the following identity: 

(50) 5 t ( I [CL + i/x / CLJ<fc) = 15e / CC**Cx d ^ 2 / C 2 C**G^x 

+tU 3 / Cldx + ^-e 3 [ C 3 CxxCxxx dx+- ae C*<Lx dx + I + J. 
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Where, 



gCxCxxxx J ^ / CxCCxxxx dx. 



One can use the Caushy-schwarz inequality to get 

(51) I / C,x \pC,x)xxx ~Z I CxxxCx 77 I CxxCxxCx 77 / ^xCxxxCx 

JS. 1 JR Z J i 1 JR 

< Mi(-|Cx| a |Cx|a + ||CE«|aK*|a + ^C^blCb), 

and 

If If 

(52) J — / CxxCCxxx o / ^xCxCxxx 

* JR 1 JR 

< Afi(||C| 2 |C BaKB |a + ^121^*10' 

for some Mi = M\(\c\ w z,<*>) > 0. If (J46J) holds, let in accordance with (J49J) M > 
be such that 

(53) \((t,x)\<M , xel, 

for as long as the solution exists. Using the Cauchy-Schwarz inequality as well as 
the fact that fi < 1, we infer from (g9|), ([50]). ifSTjl , ([52]) and ([53]) that there exists 
C(M , Mi, M, £, /u) such that 

d t £(i) < C(M , Mi, M, e, m) E(t), 

where 

= ^IC 2 + ^ K 2 + CL + ^ KLJ <*»;■ 

(Note that we do not give any details for the components of ([50)1 other than / and 
J because these components do not involve any topography term and can therefore 
be handled exactly as in [7]). An application of Gronwall's inequality enables us to 
conclude. □ 

Our next aim is to show as in the case of flat bottoms that there are solutions to 
(|24[) that blow-up in finite time as surging breakers, that is, following the pattern 
given in Proposition 3J 

Proposition 5. Let b € H°°(M.). If the initial wave profile Co G H 3 (M.) satisfies 
sup{ C o(x)}| 2 > 2 l C ^-^ + \eCl /2 ^ + \e*Cl^ 

+ 1 Co M- 1/2 + I C 1/2 Ci ^ e- 1 + | Cl' 2 Ci M - 3 / 4 s~\ 

Co= [ [Co + i^M(Co) 2 ] da > 0, (7i = \c\ w ^ > 
then wave breaking occurs for the solution of in finite time T = O(^). 

Proof. One can adapt the proof of this Proposition in the same way of the proof of 
the Proposition 6 in [7], and we omit the proof here. □ 



xei 



where 
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3.3. Well-posedness for the variable bottom KDV equation. We prove now 
the well-posedness of the equations ([29]) . ([30]) . We consider the following general 
class of equations 

, . u t + cu x + kc x u + eguu x + \[ic b u xxx = 

u u=a = u°, 

where k G R, g = § for (|30[) and g = ^ for ([29]) . This class of equations is not 
included in the family of equations stated in §3.1 because m = (here there is not 
d 2 d t u term). 

Theorem 2. Let s > | and 6 G Le£ a/so p' 6e a family of parameters 

8 = (e,/3, a,/i) satisfying \2T\ ). Assume morever that 



3c >0, V0 G p', c(ar) = y'l -/36( Q )(a;) > c . 

T/ien /or aZ/ it € i? s+1 (R), t/iere exists T > and a unique family of solutions 
(u e ) eep to bounded in C([0, f ]; i? s+1 (R)) n ^([O, f ]; H S (M)). 

Proof. As in the proof of Theorem [T] for all u smooth enough, let us define the 
"linearized" operator £(v, d) as: 

£(v,d) = d t + cd x + kc x + egvd x + \nc b dl. 

6 

We define now the initial value problem as: 
(55) (£(v,d)u = ef, 

Equation (|55p is a linear equation which can be solved in any interval of time in 
which the coefficients arc defined. We establish some precise energy estimates on 
the solution. First remark that when m = 0, the energy norm | • |x s defined in the 
proof of Theorem [T] does not allow to control the term gu xxx (for instance). Indeed, 
Lemma [5] requires m ^ to be true. We show that a control of this term is however 
possible if we use an adequate weight function to defined the energy and use the 
dispersive properties of the equation. More precisely, inspired by [5] let us define 
the "energy" norm for all s > as: 

E s {uf = \wA s u\ 2 2 

where the weight function w will be determined later. For the moment, we just 
require that there exists two positive numbers w\ , W2 such that V x G K 

W\ < w(x) < W2, 

so that E s (u) is uniformly equivalent to the standard ff s -norm. Differentiating 



e eXt E s (u) with respect to time, one gets using ([55 



2 e eXt d t {e~ EXt E s (u) 2 ) = -—E s (u) 2 - ([A s , c]d x u, wA s u) - (cdA s u,wA s u) 



-k(A s (ud x c) , wA s u) - e([A s , gv]d x u, wA s u) - e(gvdA s u, wA s u) 

—n([A s ,c 5 ]d^u,wA s u) - -n(c 5 d 3 A s u,wA s u) + e(A s f,wA s u). 
6 6 

It is clear, by a simple integration by parts that 

(56) \s(gvdA s u, wA s u) \ < eC{\v\ w i^ , \w\ W i,~ )E"(uf. 
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Let us now focus on the seventh and eighth terms of the right hand side of the 
previous identity. In order to get an adequate control of the seventh term, we have 
to write explicitly the commutator [A s , c 5 ]: 

[A s , c 5 ]d 3 x u = {A s ,c 5 } 2 dlu + Qifigu, 

where {•, -}2 stands for the second order Poisson brackets, 

{A s ,c 5 } 2 = -sd x (c 5 )A s - 2 d x + ±[sd 2 x (c 5 )A s - 2 s(s - 2)d 2 x (c 5 )A s ^d 2 ] 

and Qi is an operator of order s— 3 that can be controled by the general commutator 
estimates of (see [23] ). We thus get 

|(Qi^u,wA s m)| < CeE s {uf. 

We now use the identity A 2 = 1 — d 2 and the fact that a/3 = O(e), to get, as in 

|([ S ^(c 5 )A s - 2 - s(s - 2)d 2 x (c 5 )A s - 4 d 2 x ]d 3 x u,wA s u)\ < eC(\w\^)E s (u) 2 . 
This leads to the expression 

l/j,([A s ,c 5 ]d 3 u,wA s u) = ^u(dc 5 A s d 2 u,wA s u) + Q 2 , 
6 6 

where IQ2I < eC(\w\wi,^)E s (u) 2 . Remarking now that 

^fj l (dc 5 A s d 2 u,wA s u) = 

- S -li(d{d{c b )w)A s du,A s u) - ^(d{c b ) Wl {A s du) 2 ). 
The control of the eighth term comes in the same way: 

-/j,(c 5 d 3 A s u,wA s u) = -^i(d 3 {c 5 w)A s u,A s u) 



similarly: 



-i/i(d 2 (wc 5 )A s ,9u,A s w) - ^(d{wc 5 )A s u,A s d 2 u) 



^i(d{wc 5 )A s u,A s d 2 u) = 

^(d 2 {c 5 w)A s du,A s u) + ^ t i(d(c 5 w),(A s du) 2 



We choice here w so that 

(57) - ^(d{c 5 )w, (A s du) 2 ) + ^(d(c 5 w), (A s du) 2 ) = 

therefore if one take w = c 5 ^ -1 ^ we get eaisly ([57)) . Finally, one has 

-n([A s ,c 5 ]d x u, wA s u) + -n(c 5 d 3 A s u,wA s u) 
6 6 

= Q2 ~ ^{d(d(c 5 )w)A s du,A s u) - ^(d 3 (c 5 w)A s u,A s u) 
-^n(d 2 {c 5 w)A s du,A s u) + ^fi(d 2 {c 5 w)A s du,A s u), 



using again the fact that a/3 = 0(e) one can deduce 

'dtie-^E^u) 2 ) < e(C(E s (v)) - \)E s {uf + 4eE s (f)E s (u). 



e eXt : 
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This inequality, together with end of the proof of Theorem [IJ easily yields the 
result. □ 

4. Rigorous justification of the variable bottom equations 

4.1. Rigorous justification of the variable bottom CH equation. We restrict 
here our attention to parameters e, (3, a and /i linked by the relations 

(58) e = 0(y/jX), Pa = 0(e), [3a = 0(p 2 ). 

These conditions are stronger than (fT4"|) , and this allows us to control the secular 
effects that prevented us from proving an £P-consistency (and a fortiori a full 
justification) for the variable bottom equations of $U The notion of H s -consistency 
is defined below: 

Definition 2. Let pi be a family of parameters 9 = (e,f3,a,n) satisfying &58\ ). 
A family (C, 6 ,u )ee.pi * s H s -consistent of order s > and on [0, with the GN 
equations IfWjj . if for all 9 G p\, (and denoting h 9 = 1 + e( 9 — f3b^ a > ), 

f C 9 + [h 9 u] x = ix 2 A 

\ u% + C + eu 9 u 9 x = ~[(h e ) 3 (u e xt + eu 9 u 9 xx - e(uif)] x + i?r% 
with (r?,r|) eepi bounded in L°°([0, |], iP(M) 2 ). 

Remark 8. Since H S (M.) is continously embedded in F°°(R) for s > 1/2, the H s - 
consistency implies the L°° -consistency when s > 1/2. 

Proposition 6. Let b G H°°(M) and p G R. Assume that 

1 TP 3 1 c 9 23 
A = p. B = v h = — v , h = — p . 

y ' 1 6 2 F 6 2^ 24 

TTien: 

• For all family pi of parameters satisfying 158\) . 

• For all s > large enough and T > 0, 

• For a// bounded family (u 9 ) 0epi G C([0, j], iP(R)) solving [TUty . 
the familly (i^ 9 ,u 9 )g^ pi with (omitting the index 9) 

(59) C : = cu + + ^c 4 u xt - ep,c A \^uu xx + ^u x ], 
is H s -consistent on [0, with the GN equations US]) . 

Proof. This is clear by using the same arguments of the proof of Proposition [T] if 
we remark that now the term c x u (responsible of the secular growth effects) is of 
order 0(/i 2 ) in L°°([0, j], H S (R)). Therefore this quantity is not required in the 
identity (Ti"§|) that defines v. In the proof of Proposition [1] there is therefore no c x u 
term and the Ol=°(/^ 2 ) terms in Proposition [T] are now of order 0([i 2 ). □ 

In Proposition [6] , we constructed a family (u 9 ,(^ 9 ) consistent with the Green- 
Naghdi equations in the sense of Definition O A consequence of the following 
theorem is a stronger result: this family provides a good approximation of the 
exact solutions (u 9 , ( e ) of the Green-Naghdi equations with same initial data in the 
sense that (u e , Q 9 ) = (u 9 , C 9 ) + 0(fi 2 t) for times 0(1/ e). 
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Theorem 3. Let b G H°°(JSL), s > and pi be a family of parameters satisfying 
t58\) with (3 = 0(e). Let also A, B, E and F be as in Proposition® Lf B < then 
there exists D > 0, P > D and T > such that for all u e E H S+P (R): 

• There is a unique family (u e , C e )ee Pl € C([0, ^};H S+P (R) x H s+p - 2 ) given 
by the resolution of Mty) with initial condition Uq and formula S59\) : 

• There is a unique family (u e ,C )eepi G C([0, — ]; H S+D (R) 2 ) solving the 
Green- Naghdi equations HI 3]) with initial condition (uo,C\ t _ )- 

Moreover, one has for all 9 G pi, 

Vi 6 [0, -], \u e - u 9 |i«( [0 ,t]xR) + \C 6 ~ C e |L- ([ o, t ]xM) < Cst fi 2 t. 

Remark 9. Lt is known (see [2]^ that the Green-Naghdi equations give, under the 
scaling with f3 = 0(e) a correct approximation of the exact solutions of the 
full water-waves equations (with a precision 0(fi 2 t) and over a time scale 0(\/e)). 
Lt follows that the unidirectional approximation discussed above approximates the 
solution of the water-waves equations with the same accuracy. 

Remark 10. We used the unidirectional equations derived on the velocity as the 
basis for the approximation justified in the Theorem One could of course use 
instead the unidirectional approximation derived on the surface elevation. 

Proof. The first point of the theorem is a direct consequence of Theorem[TJ Thanks 
to Proposition^ we know that (u 6 , (f)e^ Pl is ff s -consistent with the Green-Naghdi 
equations (|13p . so that the second point of the theorem and the error estimate follow 
at once from the well-poscdness and stability of the Green-Naghdi equations under 
the present scaling (see Theorem 4.10 of [3]). □ 

4.2. Rigorous justification of the variable bottom KDV-top equation. In 

this subsection the parameters e, j3, a and fi are assumed to satisfy 

(60) e = 0(n), [3a = 0{e 2 ). 

We give first a Proposition regarding the iP-consistency result for the KDV-top 
equation. 

3 1 1 

(61) Ct + <x + + -^c b C, xxx + -c x ( = 0, 

2c D 2 

Proposition 7. Let b e H°°(R). Then: 

• For all family p\ of parameters satisfying \6th) . 

• For all s > large enough and T > 0, 

• For all bounded family {C 8 )ee P ' 1 & C([0, ^],ff s (IR)) solving 
the familly (£ e ,u e )g^ P ' with (omitting the index 9) 




is H s -consistent on [0,-^] with the equations HIS]) . 



Remark 11. Similarly, one can prove that the solution of the following KDV-top 
equations 

3 1 3 

(63) u t + cu x + -euu x + ~nc 5 u xxx + -c x u = 0, 
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071 the velocity, with 

(64) ( := cu + -u 2 - ^c 5 u xx 

4 6 

is H s -consistent with the equations Green- Naghdi equations 

Proof. One can adapt the proof of Proposition [3] in the same way as we adapted 
the proof of Proposition [T] to establish Proposition [6] (we also use the fact that if 
a family is consistent with the Boussinesq equations (|28p , it is also consistent with 
the Green-Naghdi equations (TTSJ) under the present scaling). □ 

The following Theorem deals the rigorous justification of the KDV variable bot- 
tom equation: 

Theorem 4. Let b G H°°(JSL), s > and p' x be a family of parameters satisfying 
U)Ojl with (3 = 0(e). If B < then there exists D > 0, P > D and T > such that 
for all C e G H S+P (R): 

• There is a unique family (C e , u 6 )^ G C([0, £]; H s+p - 2 (M) x H s+P ) given 
by the resolution of i61)) with initial condition £q and formula H62\ ): 

• There is a unique family (£ ,u e )0£ P > 1 G C([0, H s+D (M.) 2 ) solving the 
Green-Naghdi equations U3]) with initial condition (Co> w | t _ )- 

Moreover, one has for all 9 G <p\, 

V* G [0, |], |/ - u e U»([o >t]X K) + IC - C 8 U~([o, t] x R ) < Cst e 2 t. 

Remark 12. We used the unidirectional equation derived on the free surface ele- 
vation as the basis for the approximation justified in the Theorem [^} One could of 
course use instead the unidirectional approximation i63\) derived on velocity. 

Proof. The first point of the theorem is a direct consequence of Theorem[5J Thanks 
to Proposition^ we know that (u e , £ e )g ep ' i is iP-consistent with the Green-Naghdi 
equations p3[) . so that the second point follows as in Theorem [3] □ 

Acknowledgments. The author is grateful to David Lannes for encouragement 
and many helpful discussions. 
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